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ABSTRACT 



The physical interpretation of the spectral line polarization produced by the 
joint action of the Hanle and Zeeman effects offers a unique opportunity to obtain 
empirical information about hidden aspects of solar and stellar magnetism. To 
this end, it is important to achieve a complete understanding of the sensitivity of 
the emergent spectral line polarization to the presence of a magnetic field. Here 
we present a detailed theoretical investigation on the role of resonance scattering 
and magnetic fields on the polarization signals of the Ball Di and D2 lines of the 
Fraunhofer spectrum, respectively at 4934 A and 4554 A. We adopt a three-level 
model of Ba II, and we take into account the hyperfine structure that is shown 
by the 135 Ba and 137 Ba isotopes. Despite of their relatively small abundance 
(18%), the contribution coming from these two isotopes is indeed fundamental 
for the interpretation of the polarization signals observed in these lines. We con- 
sider an optically thin slab model, through which we can investigate in a rigorous 
way the essential physical mechanisms involved (resonance polarization, Zeeman, 
Paschen-Back and Hanle effects), avoiding complications due to radiative trans- 
fer effects. We assume the slab to be illuminated from below by the photospheric 
solar continuum radiation field, and we investigate the radiation scattered at 
90°, both in the absence and in the presence of magnetic fields, deterministic 
and microturbulent. We show in particular the existence of a differential mag- 
netic sensitivity of the three-peak Q/I profile that is observed in the D2 line in 
quiet regions close to the solar limb, which is of great interest for magnetic field 
diagnostics. 
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1. INTRODUCTION 

Probably, the most interesting aspect of spectropolarimetry is that it allows us to diag- 
nose magnetic fields in astrophysics. To this end, it is crucial to achieve a complete physical 
understanding of the magnetic sensitivity of the emergent spectral line radiation given the 
fact that it can occur through a variety of rather unfamiliar physical mechanisms, not only 
via the Zeeman effect (e.g., Landi Degl'Innocenti & Landolfi 2004 (hereafter LL04); see also 
the reviews by Stenflo 2003 and Trujillo Bueno 2003). In this respect, the main aim of this 
paper is to help decipher the physical mechanisms that control the magnetic sensitivity of 
the polarization of the D-lines of Ba II, with particular interest in developing a powerful 
diagnostic tool for mapping the magnetic field of the lower solar chromosphere. 

In the atmospheres of the Sun and of other stars there is a fundamental mechanism 
producing polarization in spectral lines, which has nothing to do with the familiar Zee- 
man effect. There, where light escapes through the stellar "surface", the atomic system is 
illuminated anisotropically. The radiative transitions produce population imbalances and 
quantum coherences between pairs of magnetic sublevels, even among those pertaining to 
different levels. The mere presence of this so-called atomic level polarization produces spec- 
tral line polarization, without the need of a magnetic field. This is usually referred to as 
resonance line polarization. The important point is that a magnetic field can modify the 
atomic polarization of the upper and/or lower levels of the spectral line under consider- 
ation and the ensuing polarization of the emergent spectral line radiation. Interestingly, 
the possible presence of crossings and repulsions among magnetic sublevels of fine-structure 
and/or hyperfine-structure multiplets can enhance dramatically the magnetic sensitivity of 
the emergent spectral line polarization. A remarkable example is the enhancement of the 
line-core scattering polarization of the D 2 line of Na I by a vertical magnetic field, which is 
due to interferenc es between particular hyper fine structure (HFS) magnetic sublevels of the 



2 P 3 / 2 upper level (ITrujillo Bueno et al.l 120021) . It is of interest to note that this theoretical 



prediction was observationally confirmed by lStenflo et al.l (120021 ) via filter polarimetry of the 
solar atmosphere. 

In contrast with the case of sodium, which has one single isotope with nuclear spin 
1=3/2, barium has five (even) isotopes with 1=0 (with an overall abundance of 82.18%) and 
two (odd) isotopes with 1=3/2 (with an abundance of 17.82%). Moreover, the HFS split- 
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ting of the odd isotopes of barium is about a factor five larger than for the case of sodium. 
Obviously, the emergent fractional linear polarization (i.e., the Q/I profile, where I and Q 
are two of the St okes paramete rs) has contributions from all the barium i sotopes. In fact, 
as poi nted out by IStenflol (119971 ). the Q/I pattern of the Ball D 2 line that IStenflo fc Keller 
( 119971 ) observed in very quiet regions close to the solar limb shows a three-peak structure, 
with a prominent central Q/I peak due to the even isotopes (which are devoid of HFS) and 
two less significant peaks in the red and blue wings caused by the contributions from the odd 
isotopes (which have HFS). Therefore, we think that for the D 2 line of Ball we should also 
have enhancement of scattering polarization by a vertical field, but only around such wing 
wavelengths because the required interferences occur only between the magnetic sublevels of 
the 2 i"3/2 upper level of the barium isotopes endowed of HFS. Actually, the scientific moti- 
vation that led us to undertake the theoretical investigation presented here was to develop 
a novel plasma diagnostic tool based on the idea that such isotopes of barium must have a 
different behavior in the presence of a magnetic field, with respect to those devoid of HFS. 

While the physical origin of the observed Q/I profile of the Ba II D 2 line seems to 
be clear, nobody h as yet been able to model the Q/I profiles of the Ball Dx line that 
Stenflo et al.l (119981 ) observed in two diffe rent regions c l ose to the solar limb. Interestingly, 
while the Q/I profile shown in Fig. 1 of IStenflo et al.l (119981 ) might perhaps be the result 
of the interaction of the Di line with the continuum (see the bottom panel to the right- 
hand-side), the symmetric Q/I profile shown in the panel just above the previous one is ver y 
similar to the enigmatic Q/I profile of the Nal Di line (see also Fig. 3 in IStenflo et al.ll2000l ). 
Although the main scientific target of our paper is to understand the magnetic sensitivity of 
the above-mentioned Ball D 2 line, we present also some results of our Q/I calculations for 
the Di line with the aim of helping to clarify its physical origin. 

The outline of this paper is the following. The formulation of the problem is presented in 
Section [21 where we establish our modeling assumptions, we briefly discuss the relevant equa- 
tions, and we describe the atomic model, showing the behaviour of the magnetic sublevels 
of the odd barium isotopes in the presence of an increasing magnetic field. The magnetic 
sensitivity of the atomic polarization of the lower and upper levels of the Ba II D-lines is 
discussed in Section |3l pointing out the similarities with the behavior of the Nal levels. Sec- 
tion H] focuses on the emergent spectral line polarization in the absence and in the presence 
of a magnetic field, with emphasis on the D 2 line of Ball, but showing also some interesting 
results for the Di line. Finally, Section [5] summarizes our main conclusions with an outlook 
to future research. The three appendices give detailed information that may help the reader 
to understand better the complexity of the problem we are investigating. 
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2. 



FORMULATION OF THE PROBLEM 



We consider an optically thin slab of Ba II ions illuminated from below by the pho- 
tospheric solar continuum radiation field (see Fig. [[]). The slab is assumed to be located 
1000 km above the r 50 oo = 1 photospheric level, the approximate height at which, according 
to semi-empirical models of the solar atmosphere, the line-core optical depth of the D 2 line 
is unity along a line-of-sight specified by /i = cos# = 0.1, with 9 the heliocentric angle. 

In the following subsections we summarize the basic equations derived within the frame- 
work of the quantum theory of spectral line polarization that will be used in this work (§ 12. ip . 
we describe the atomic model that we have adopted to describe the Ball ion (§ 12.21) . and fi- 
nally we show how the main properties of the continuum photospheric radiation field incident 
on the slab have been calculated (§ 12.31) . 



We describe the excitation state of the Ball levels using the density matrix formalism, 
a robust theoretical framework very suitable to treat the atomic polarization (population 
imbalances and quantum interferences, or coherences, among the magnetic sublevels) that 
anisotropic pumping processes induce in an atomic system. Referring to one of the isotopes of 
barium (see § 12.21) that shows HFS, we indicate with / its nuclear spin quantum number. In 
the absence of magnetic fields, using Dirac's notation, the energy eigenvectors can be written 
in the form \a J I F f >, where a represents a set of inner quantum numbers (specifying the 
electronic configuration and, if the atomic system is described in the L-S coupling scheme, 
the total electronic orbital and spin angular momenta), J is the total electronic angular 
momentum quantum number, F is the total angular momentum quantum number (electronic 
plus nuclear: F=J+I), and / is its projection along the quantization axis. In principle, 
to have a suitable description of atomic polarization, one has to take into account all the 
coherences of the form 



where p is the density matrix operator. The approximation that is used in this paper consists 
in restricting the description of the atomic system to the J-diagonal density matrix elements 



or, in other words, in neglecting coherences between different J-levels. The resulting model- 
atom is referred to as the "multi-level atom with hyperfine structure" according to LL04. 
The approximation appears to be fully justified for the investigation of the Ba II D-lines 
given the large frequency separation between the levels 6p 2 P±/2 and 6p 2 P3/2 (see § 12 . 2[) . and 



2.1. The Basic Equations 



< a J I F f \ p\a J' I F' f > , 




< a J I F f \ p\a J I F' f > , 



(2) 
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Fig. 1. — Geometry of the problem under investigation. 

given the relatively low abundance of barium in the solar atmospher^]. 

In the presence of a magnetic field, according to the general approach of the Paschen- 
Back effect theory, the energy eigenvalues and eigenvectors have to be found by diagonaliza- 
tion of the total Hamiltonian (unperturbed atomic Hamiltonian plus magnetic Hamiltonian) 
on each J-level subspace. Taking the quantization axis along the magnetic field direction, 
it can be demonstrated that the total Hamiltonian commutes with the projection along the 
quantization axis of the total angular momentum operator (/ is a good quantum number), 
while, in general, it does not commute with the total angular momentum operator. The 
eigenvectors of the total Hamiltonian can be expressed in the form (e.g. LL04) 



where the index i labels the energy eigenstates belonging to the subspace corresponding to 
assigned values of the quantum numbers a, J, I and /, and where the coefficients C l F can be 
chosen to be real. In the energy eigenvectors representation the atomic system will therefore 
be described by means of the matrix elements 



If the magnetic field is so weak that the magnetic energy is much smaller than the energy 
intervals between the HFS F-levels, we are in the so-called Zeeman effect regime (of HFS), 



1 Note that this approximation is not justified for the D-lines of Nai and for the H and K lines of Can, 
given the larger abundance of these ions. 




(3) 



F 




(4) 
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where the energy eigenvectors are still of the form \a J I F f > (C F ~ Spi), and the splitting 
between the HFS magnetic sublevels is linear with the magnetic field strength. For stronger 
magnetic fields it is necessary to apply the Paschen-Back effect theory, and one enters the 
so-called incomplete Paschen-Back effect regime. In this regime the energy eigenvectors have 
the general form of equation ([3]) (the magnetic field produces a F-mixing of the various HFS 
levels originating from a particular J-level) and, as we will show in detail for the case of Ball 
in § \2.2\ the splitting among the various HFS magnetic sublevels is no longer linear with the 
magnetic field. Several crossings among HFS magnetic sublevels with different / quantum 
number take place in this regime, as well as a repulsion among the magnetic sublevels 
with the same / quantum number. This behaviour of the magnetic s ubley els has important 



consequences on the atomic polarization, as pointed out by iBommierl (1198CH ) and described in 
detail in LL04, and produce interesting effects, sometimes referred to as level crossing effect 
and anti-lev el- cms sing effect, on the polarization signals produced by resonance scattering 



(e.g. 



Trujillo Bueno et al.ll2002l ). In § 14.31 we will show their effect on the linear polarization 
of the Ba II D2 line. If the magnetic field strength is further increased the so-called complete 
Paschen-Back effect regime is reached. In this regime the energy eigenvectors are of the form 
\a J I Mj Mj >, and the splitting among the HFS magnetic sublevels is again linear with the 
magnetic field strength: the atom behaves in this regime as if it were devoid of HFS. Going 
from the Zeeman effect regime to the complete Paschen-Back effect regime, the magnetic 
field produces therefore an energy eigenvectors basis transformation. 

In the following we will work in the spherical statistical tensor representation. The 
conversion of the density matrix elements of equation (J3J) into this representation is given by 
the relation (cf. LL04) 

PaJi(if, i'f) = ^Q(aJ/, /) C F \{a J I, f)p aJI (Ff, F' f ) = 

FF' 

= j2cF^ji,f)cu^ji,f)J2(- 1 ^ fV ^^( F f S K ) aJIp Q^ F,) • (5) 

FF' KQ V J J V/ 

The Statistical Equilibrium Equations (SEEs) and the radiative transfer coefficients for a 
multi- level atom with HFS, in the spherical statistical tensor representation, written taking 
the quantization axis directed along the magnetic field, can be found in § 7.9 of LL04. Here 
we write only the expression for the emission coefficient in the transition between the upper 
level (a u , J u ) and the lower level (ag, Je) 

2hu 3 hv 

£j (u, O) = ——N(2J U + l)B(a u JJ -> a e J t I) 
c z Air 

x V / 3(2F + 1)(2F U +1) E (~l) 1+K ~ fuW 

KQK U Q U i u F u F' u F"ii Ff Fi fufUill 1 
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(6) 



where j = 0, 1, 2, 3, standing respectively for the Stokes parameters /, Q, U and V, Af is 
the number density of atoms, B(a u J u I — > a^J^I) is the Einstein coefficient for stimulated 
emission, Tqiji fi) is a geometrical tensor (cf. LL04), and $ the profile of the line. 

It is important to note that the previous equations are valid under the flat-spectrum 
approximation. For a multi-level atom with HFS this approximation requires that the inci- 
dent radiation field should be flat (i.e. independent of frequency) across a spectral interval 
Av larger than the frequency intervals among the HFS levels (possibly split by the magnetic 
field), and larger than the inverse lifetimes of the same levels. This is a good approximation 
for the Dx and D 2 lines of Ba II if we restrict to magnetic fields smaller or of the order of 
1 kG. 



2.2. The Atomic Model 

We adopt a three- level model of Ball consisting in the ground level (6s 2 S1/2), the upper 
level of the Dx line (6p 2 P\/2) and the upper level of the D 2 line (6p 2 F 3 / 2 ). There are seven 
stable isotopes of barium, whose mass numbers and relative abundances are listed in Table 1. 
In this work we take into account the contributions coming from all the seven isotopes. 

The mass and volume differences between the nuclei of the various isotopes involve small 
but appreciable differences on the energies of the fine structure levels of different isotopes 
(isotopic effect). We use the values of the isotopic shifts in the D x and D 2 lines listed in 
Table 1 to correct the energies of the 2 Fi/ 2 and 2 F 3 / 2 leve ls of the variu os isotopes. For the 



reference isotope (138) we use the energy values given by iMoord ( 119581 ). 

Isotopes with even mass number have nuclear spin 1 = 0, while those with odd mass 
number (135 and 137) have nuclear spin 1 = 3/2. The odd isotopes show therefore HFS due 
to nuclear spin. Introducing the total angular momentum, characterized by the quantum 
number F, we observe (see Fig. |2]) that the levels 2 5*i/ 2 and 2 Fi/ 2 split into two HFS levels 
(F = l,2), while the level 2 F 3/2 splits into four HFS levels (F = 0,1,2,3). It is possible to 
demonstrate that the HFS Hamiltonian can be expressed as an infinite series of electric and 
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Isotope 


Abund. 

(%) a 


I 


Isotope Shifts (MHz) 1 


HFS Constants (MHz) 2 


Di 


D 2 


2 S\/2 2 Pl/2 

A A 


A B 


130 Ba 


0.106 





355.3 b 


372. 3 b 






132 Ba 


0.101 





278.9 b 


294.9 b 






134 Ba 


2.417 





222.6 C 


234.6 C 






135 Ba 


6.592 


3/2 


348.6 b 


360.7 b 


3591.67 d 664.6° 


113.0° 59.0° 


136 Ba 


7.854 





179.4 b 


186. 9 b 






137 Ba 


11.232 


3/2 


271. l b 


279.0 b 


4018.87 d 743.7° 


127.2° 92.5° 


138 Ba 


71.698 





reference isot. 







1 A positive I.S. means that the line is shifted to higher frequencies with respect to the reference isotope. 

2 The HFS constant B is defined accord i ng to t he convention of the A merican literature. 

a NIST on-line database; IWendt et ail |l984h : iWendt et all |l988h : ^Becker et ail fjl98lh : ^Villemoes etal 

(|l993h . 



Table 1: Isotopes considered in this work. 



magnetic multipoles (e.g. iKopfermannl Il958l ) . To calculate the energies of the various HFS 
levels we consider the first two terms of the series (magnetic dipole and electric quadrupole 
terms), and we use the HFS constants listed in Table 1. This investigation will clearly show 
the importance of the HFS effects for a correct modeling of the polarization produced by 
scattering processes in a stellar atmosphere. 

The Lande factors have been calculated theoretically, assuming L-S coupling for the 
Ball ion. The values obtained differ by less than 1.3% from the experimental ones ((Moore 
19581 ). reported in Figure EJ This can be considered as a proof that the L-S coupling is quite 
a good approximation for the Ball ion. 

The energies of the HFS magnetic sublevels of the isotope 137, as functions of the 
magnetic field strength, are shown in Figure [2j We can see that in the range between 
and 1000 G the splitting of the various magnetic sublevels originating from the ground level 
2 S\/2 is linear with the magnetic field strength (Zeeman effect regime). A similar behaviour is 
shown by the magnetic sublevels originating from the level 2 P\/2 for magnetic fields smaller 
than about 600 G. For stronger magnetic fields the linearity of the splitting appears to 
be slowly lost, which indicates that the incomplete Paschen-Back effect regime is reached. 
The splitting observed among the magnetic sublevels originating from the level 2 P^/2 shows 
instead that a complete transition from the Zeeman effect regime to the complete Paschen- 
Back effect regime takes place for magnetic fields ranging from to 1000 G. As described 
in 



12.11 in the intermediate incomplete Paschen-Back effect regime several level crossings 
among HFS magnetic sublevels can be observed, as well as a repulsion among the sublevels 
with the same / quantum number. 



- 9- 



Energy (MHz) 



658117334.0 



607426618.6 



g=1.32 



g=0.672 



g= 1.974 



ISOTOPE 137 

F Energy (MHz) 

— 3 658117643.4 

— 2 658117169.3 



658117007.4 
658116972.7 



607427176.4 



607425689.0 



3014.2 



-5023.6 





1000 


N 


500 


X 









Ld 


-500 




-1000 




1000 


N 


500 


X 


5^ 





Ld 


-500 




-1000 




200 400 600 800 1000 



200 400 600 800 1000 





4000 






2000 






N 








X 











-2000 








-4000 








-6000 





200 400 600 800 1000 
B (G) 



Fig. 2. — Energies, corrected for the isotopic effect, of the fine-structure and HFS levels of the isotope 
137 Ba II. The right panels show the energies of the HFS magnetic sublevels as functions of the magnetic 
field strength (in each panel the zero of the energy scale is chosen at the energy of the corresponding fine 
structure J-level). 



2.3. The Incident Radiation Field 

As already stated in § [2J, we consider an optically thin plane-parallel slab, composed of 
Ba II ions, located at approximately 1000 km above the 75000 = 1 photospheric level, and we 
assume that the slab is illuminated from below (hence, anisotropically) by the photospheric 
continuum radiation. Under these hypothesis the atomic polarization can be calculated solv- 
ing directly the SEEs for the given continuum radiation field coming from the photosphere. 

Let us take a reference system with the z axis directed along the local vertical, and let 
us describe the continuum radiation field incident on the slab by means of the tensor 

3 

J a— n 
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Line 


A (A) 




n v 




Di 


4934.09 


0.955x10 s 


0.323xl0~ 2 


0.159 


D 2 


4554.03 


1.17xl0 8 


0.225xl0~ 2 


0.176 



Table 2: Wavelength (in air), and Einstein coefficient of the transitions considered; mean number of photons 
and anisotropy factor of the photospheric continuum at the wavelength of the transitions, 1000 km above 
the T5000 = 1 level. 



where Si = I, Q, U, V. Assuming that the incident radiation field is unpolarized and has 
cylindrical symmetry around the local vertical, it is easy to verify that the only non vanishing 
components are 

J °(,)^/(^) and J 2 (z,)=y^(^=(3/i 2 -l)/(^ M )) . (8) 

Note that J° is just the mean intensity of the incident radiation (averaged over all directions), 
while Jq gives a measure of the anisotropy of the radiation fielcjE Instead of Jq and J5 2 , we 
can use the quantities n(v), the mean number of photons, and w u , the so called anisotropy 
factor. The new quantities are related to the previous radiation field tensor components 
through the relations 

»M = ^oV> and «V = ^|g- (9) 

To calculate the values of n(y) and w v of the photospheric continuum at the height of 
1000 km above the visible solar "surface", at the frequencies of the Di and D2 lines of Ball, 
we follow § 12.3 of LL04. The values of the specific intensity of the rad iation coming from 



the solar disk center and of the limb- darkening coefficients are taken from lPiercd (120001 ) . The 
values obtained for n u and w u are listed in Table 2. At this point the SEEs can be solved 
numerically. Their expression becomes simpler if we rewrite them in the reference system 
with the quantization axis directed along the local vertical direction, as in this case only two 
components of the radiation field tensor are non zero (Jq and Jq). It can be demonstrated 
that all the radiative rates are invariant under a rotation of the reference system so that only 
the magnetic kernel has to be modified with respect to the expression given in equation (7.66) 
of LL04. 



2 In particular Jq quantifies the unbalance between vertical and horizontal illuminations. 
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3. THE POLARIZATION OF THE ATOMIC LEVELS 

We solve numerically the SEEs (which implies, for each isotope with HFS of our model 
atom, the solution of a linear system of 384 equations in the unknowns Pq{F, F')) for mag- 
netic field strengths between and 1000 G, and for various inclinations of the magnetic 
field with respect to the local vertical. We recall that the p®(F,F) elements quantify the 
populations of the various F-levels, the Pq elements {alignment components) contribute to 
the linear polarization of the scattered radiation, while the pg elements (orientation compo- 
nents) contribute to the circular polarization of the scattered radiation. As shown in LL04, 
an anisotropic, unpolarized, flat spectrum radiation field generally induces only alignment 
in the atomic system, while orientation can be originated by the so called alignment-to- 
orientation conversion mechanism. Note that all the levels of the isotopes 135 and 137, 
because of the HFS, can carry alignment while, for all the other isotopes, only the level 
2 P 3 / 2 , the upper level of the D 2 line, can carry it. 

Let us consider isotope 137. In complete analogy wi th the case of Nal , having a single 



isotope with 7 = 3/2 (see Trujillo Bueno et al. 2002 and lCasini et al.ll2002L and the discus- 
sion therein), only the level 2 P^/2 is polarized directly via the anisotropic illumination. The 
ground level becomes polarized because of a transfer of polarization via spontaneous emission 
in the D 2 line, while the level 2 P\/2 becomes polarized via radiative absorption (repopulation 
pumping) in the Di line. This explains the fact that the upper and lower levels of the Di 
line are equally sensitive to the magnetic field strength, independently of its inclination (see 

Fig. ED. 

It is well known that a magnetic field is able to modify the atomic polarization, and 
therefore the polarization of the scattered radiation (Hanle effect^]). The behaviour of the 
various spherical statistical tensors, written in the local vertical reference system, as func- 
tions of the magnet ic field inclination and str ength, is qualitatively equal to the case of 



Nal, investigated by lTrujillo Bueno et al.l (120021 ). We point out that the first decrease of the 



atomic polarization of the Di levels, which takes place only for inclined fields (see Fig. [3]), is 



3 In this work by Hanle effect we mean any modification of the atomic polarization which is due to the 
action of a magnetic field. Note, however, that in the literature it is often meant by Hanle effect only the 
relaxation of coherences (defined in the magnetic field reference system) having Q ^ 0. Within this second 
meaning, it is often stated that there is no Hanle effect in the presence of a vertical magnetic field. Note 
that according to our definition this statement is true only if we are dealing with an isolated level (i.e. if we 
neglect the quantum interferences among the magnetic sublevels originating from different hyperfine (or fine) 
structure levels). We prefer to adopt the former more general definition because, as we are dealing with a 
quite complex atomic system with HFS, and as we are investigating the role of magnetic fields with complex 
configurations (random-azimuth and microturbulent), it becomes quite difficult to understand which effects 
due to the magnetic field can be considered as 'Hanle effect' according to the latter definition. 
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Fig. 3. — Fractional atomic alignment, cr§ = pg(F, F) / p®(F, F) , of the various HFS levels of the isotope 
137, calculated in the local vertical reference system, as functions of the magnetic field vector, for various 
inclination angles. Top row 6>s=0°, middle rows #b = 30° and #b = 60°, bottom row #b = 90°. The angle 9b 
is defined in Figure [TJ 



due to the so called lower-level Hanle effecto As can be demonstrated (see LL04 for details), 
the spherical statistical tensors Pq{F, F) (that describe population imbalances and quantum 
interferences among magnetic sublevels originating from the same F-level, and that mainly 
affect the polarization of the scattered radiation at the line center) are significantly modified 
by a magnetic field when the Zeeman splitting is of the same order of magnitude as the 
inverse lifetime of the level. That is, as a rough estimation, when the magnetic field ranges 
between (see LL04) 

0.1 B c < B < 10 B c , (10) 



4 This definition has been formulated within the latter definition of the Hanle effect given in the footnote[3] 



13 



with 

me 9h 

where tuf e (in seconds) is the radiative lifetime of the lower or upper level of the line transi- 
tion under consideration, is its Lande factor and B c is the critical magnetic field intensity 
in Since the relevant atomic level here is the ground level, it is important to note that 
its radiative lifetime is tm e ~ l/(Be u Jq). As seen in Figure [3] this first decrease takes place 
for magnetic fields of the order of 10 _1 G, consistently with our simpli fying as s umpt ion that 



the pumping radiation field is the continuum radiation tabulated by IPiercd (120001 ) . If, on 



the contrary, one takes into account the line profile (resulting in a sm aller value of J°) the 



decre ase will take place for smaller magnetic fields (cf. the Nal results of lTrujillo Bueno et al. 



20021 ). Note that this decrease of the atomic polarization in the ground level has a feedback 
even on the D2 upper level. Concerning the second sudden decrease of the atomic polariza- 
tion of the Di levels, for B larger than 100 G, we rec all that this is due to the in hibition 



of the repopulation pumping mechanism discussed by iTrujillo Bueno et al.l (120021 ) and by 



Casini et al.l (120021 ). which sets in when the electronic and nuclear angular momenta, J and 
/, of the 2 Ps/2 level are decoupled, (i.e. when this level, the only one that can carry alignment 
even in the absence of HFS, enters the complete Paschen-Back effect regime). Note that for 
the case of sodium this sudden decrease occurs for B larger than 10 G simply because the 
complete Paschen-Back effect regime is reached for weaker magnetic fields in sodium than 
in barium. 

As expected, because of the symmetry of the problem, for a vertical magnetic field only 
the components with Q = are non zero. For different orientations of the magnetic field we 
have in general contributions coming from all the density matrix elements; in particular it 
is possible to demonstrate that the components with Q = are independent of the magnetic 
field azimuth, the components with Q = 1 change sign under an azimuth rotation of 180°, 
the components with Q = 2 change sign under an azimuth rotation of 90°, and so on. 

The expressions of the emission coefficients (eq. [H]), as well as the expressions of all the 
other radiative transfer coefficients given in LL04 hold in the magnetic reference system. 
Therefore we have to transform the spherical statistical tensors, obtained solving the SEEs 
written in the local vertical reference system, into the magnetic field reference system. In- 
dicating with [pq(F,F')]b the spherical statistical tensor components in the magnetic field 
reference system, and with [pq(F,F')]v the spherical statistical tensor components in the 
local vertical reference system, we have 

lp%(F,Fi)] B = J2ipU F > F ')]v D %>Q( R r > (12) 



5 The previous expression of B c is exact only for an isolated level (see footnot£[3]) . 
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where Dq,q(R) is the rotation matrix calculated for the rotation R which carries the ver- 
tical reference system into the magnetic reference system (referring to Fig. [I] we have 
R — (xb,#b,0)), and where the apex "*" indicates the complex conjugate. We observe 
that after the rotation, because of the symmetry of the problem, the spherical statistical 
tensors in the magnetic reference system do not depend on the azimuth but only on the 
inclination of the magnetic field with respect to the local vertical reference system. 



4. THE POLARIZATION OF THE EMERGENT SPECTRAL LINE 

RADIATION 

For the case of a tangential observation in a plane-parallel atmosphere, it can be shown 
that, under the approximation of a weakly polarizing atmosphe re {ej ^> Eq, en, ey; m 3> 



VQjVu,Vv), the emergent fractional polarization is given by (e.g.. iTrujillo Buenol 120031 ) 
X(v,ti) e x {v,Q) r} X {y,to) 



I(v,Q) ei(u,Q) 771(1/, n) 



with X = Q,U,V . (13) 



We point out that the first term of equation ([TBI is the contribution to the emergent fractional 
polarization due to selective emission processes, while the second one is caused by dichroism 
(selective absorption of polarization components). As shown in Figure [31 for the D 2 line the 
contribution due to dichroism is much smaller than that due to selective emission. For this 
reason, from now on we will describe the polarization properties of the radiation emergent 
from the slab using the relation 

X(v,fl) _ e x (u,Cl) 

i(y,n) £,(!/, n) ' 1 ) 

We recall that the polarization properties of the emergent radiation will always be described 
assuming the reference direction for positive Q parallel to the slab. 

It is important to remember that the expressions for the emission and absorption co- 
efficients given in LL04 take into account only the line processes. For this reason, in order 
to be able to reproduce qualitatively the observed profiles, we need to add the contribution 
coming from the continuum. Assuming that the continuum is not polarized and constant 
across the line, we have 

1(1/, n) ej(i/,n) + e / ' 1 1 

where ef is the continuum contribution to the intensity of the emergent radiation, and where 
the apex " I " is to recall that the corresponding quantity refers only to the line processes. 
In Appendix [S] we show and briefly discuss some results obtained by applying equation (115]) . 
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The Ba II D-lines that we are investigating are strong lines: according to theoretical 
models of the solar atmosphere the wings of these lines originate in the photosphere, while 
the line-cores in the high photosphere-low chromosphere. The optically thin slab model il- 
luminated by the photospheric continuum that we are considering in this paper is therefore 
just a zero-order approximation. Nevertheless, it allows us to take into account in a very 
rigorous way the atomic physics involved in the problem, and to understand its essential role 
on the magnetic sensitivity of the polarization profiles of these lines, avoiding complications 
coming from radiation transfer effects. This is the first step of our investigation, in a forth- 
coming paper we will propose more realistic models, where radiation transfer effects will be 
taken into account. 

Once the SEEs in the vertical reference system have been solved numerically, and the 
spherical statistical tensors have been rotated to the magnetic field reference system, we can 
calculate the emission coefficients for 90° scattering (9 = 90° and x = 0° i n Fig- [I]) by means 
of equation (jSJ), and the polarization of the scattered radiation through equation (fT5l) . In 
the following subsections we present our results for the Ba II D 2 and Di lines, for various 
magnetic field configurations. 



4.1. D 2 Line - No Magnetic Field Case: 
Origin of the Three Peaks Structure and Choice of Parameters Values 

As a first step, we observe the laboratory positions of the various HFS components 
of the D 2 line. The three isotopes without HFS contribute to the D 2 line with just one 
component each, those which have HFS contribute to this line with six components each, 
which may overlap. As seen in Figure H] (right panel), it is possible to divide the various 
HFS components into three groups. The central group, at about 4554. 03A, is composed by 
the five components (three visible in the figure) due to the five isotopes without HFS (note 
that the main contribution comes from the isotope 138 because of its high abundance, while 
isotopes 130 and 132 bring a negligible contribution, not visible in the figure). The other 
components, due to the isotopes with HFS, fall at different wavelengths but can be gathered 
into two groups because of the large splitting of the ground level into the two F = 1 and 
F = 2 HFS levels (see Fig. [2]). In particular, the group at about 4553. 995A is composed by 
the HFS components, of both isotopes 135 and 137, associated to the transitions towards 
the lower F = 1 HFS level, while the group at about 4554. 045A is composed by the HFS 
components associated to the transitions to the lo wer F = 2 HFS level. As we will see in 



more detail below, and as already pointed out by IStenflol (119971 ) , the origin of the three 
peaks structure of the D 2 line lies in this splitting of the various components into these three 
groups. Similar considerations could be done about the position and relative strength of the 
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Fig. 4. — Laboratory positions of the various HFS components of the Di (left panel) and D2 (right panel) 
lines resulting from the seven isotopes considered. The horizontal axis gives the wavelength (in A) measured 
from 4934.075 A for the Di line, and from 4554.029 A for the D2 line. These are the wavelengths (in air) of 
the isotope 138 Di and D 2 lines, respectively. 

various components of the Di line, with the difference that now the upper level just has two 
HFS levels (instead of four). 

Let us begin our analysis of the D 2 profile by considering only the isotope 138 (without 
HFS), and let us assume a Doppler width of 30 mA. For this isotope, as shown in the left 
panels of Figure [51 the ratio £n{\) / e\(X) (often referred to as fractional polarization) is con- 
stant and different from zerqj. Adding the contribution of the continuum, the same ratio 
remains unchanged in the line-core (where e \ ^> ef) , while it goes to zero at the wavelengths 
corresponding to the wings of the intensity profile. The same considerations can be done for 
all the other isotopes without HFS. Let us consider now the isotope 137 (with HFS). For 
this isotope, assuming the same Doppler width, the profiles of the line emission coefficients 
ej and s 1 q show two peaks at the wavelength positions of the two groups of HFS transitions 
(see Fig. [5j). The ratio £q(A)/£j(A) is no longer constant but decreases showing two broad 
minima at the wavelength positions corresponding to the wings of the e\ and Eq profiles, and 
assumes the same value as the isotopes without HFS moving away from the line-core^. This 



6 See Appendix 151 for an analytical proof of this result. 
7 See Appendix iBl for an analytical proof of this result. 
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Fig. 5. — Theoretical profiles of the emission coemcients ej(X) and Sq(X), and of the ratio Eq{\) / {e\{\) + 
ef), of the Ba n D2 line, for the isotope 138, for the isotope 137 and for all the seven isotopes together, in 
the absence of a magnetic field. The emission profiles are normalized to the maximum value of the intensity 
emission profile, calculated taking into account the contribution of all the seven isotopes, (£j) max - The 
vertical dashed lines show the positions of the three groups of transitions (see text). The last row shows 
the Eq{\) / {e\{\) + ef) profiles, without continuum (solid), with a continuum ef/(e l I ) niax of 10 -5 (dot), of 
3 x 10" 5 (long dash), and of 9 x 1(T 5 (dash-dot). 



profile clearly shows the depolarizing effect of the HFS. The role of the continuum is the 
same as observed for the isotopes without HFS. The same arguments hold for the isotope 
135. 

This investigation on the isotopes 138 and 137 shows that the central peak of the ob- 
served Q/I profile is due to the isotopes without HFS, while the two secondary peaks are due 
to the isotopes with HFS. The position and amplitude of these secondary peaks appear to 
be strongly dependent on the background continuum emissivity ef, a physical quantity that, 
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given the exploratory character of this paper, we will just parametrize in order to reproduce 
at best the observed profile. 

We consider now all the isotopes together and we adjust the Doppler width, the anisotropy 
factor an d the continuum intens ity in order to obtain the best fit to the Q/I profile ob- 
served by lStenflo fc Keller! (119971 ). still assuming that no magnetic field is present. Changing 
the Doppler width we can modify the separation between the two minima of the profile 
e£(A)/ej(A). To obtain the same separation as the observed profile we need a value of about 
30 niA. This value seems to be very reasonable as it can be obtained assuming a temperature 
of about 6000 K and a microturbulent velocity of about 1.8 km/s, values which are in good 
agreement with those given by semi-empirical chromospheric models at the height of about 
1000 km. 

Modifying the anisotropy factor we simply scale the ratio eh(X) / e\(\) . Radiative trans- 
fer effects are disregarded in our model and, since the Ba II D 2 line is a strong line, we 
can expect that the calculated value for the anisotropy factor taking into account only the 
photospheric continuum (see § 12.31) . is probably overestimated. We find indeed that the 
anisotropy factor has to be decreased to the value of 0.037 (approximately 1/5 of the value 
0.176 mentioned in § 12. 3p to obtain the observed value of the ratio Q/I at the wavelength 
position of the central peak. 

In this paper we are not taking into account the continuum processes in a rigorous way. 
As mentioned above, we describe the effect of the continuum simply through the parameter 
Ej. As shown in Figure the continuum modifies the wings of the profile. Taking as a ref- 
erence the "red" secondary peak of the observed profile, we find that the best fit is obtained 
assuming a value for ej/(e l T ) maiX of 9 x 10 -5 . With these values of the parameters, applying 
equation (fl5|) . we get the profile shown in the panel b of Figure [6] which reproduces quite 
we ll the observed profile, and which is very similar to a theoretical profile already obtained 
bv lstenflol Jl997h . 

The main aim of this work is to investigate the magnetic sensitivity of the linear and 
circular polarization of the D-lines of Ball, and not merely to reproduce as better as possible 
the observed profiles of these lines. For this reason, we prefer to perform our investigation by 
sticking to the values of the parameters found for the simpler, unmagnetized case. Obviously 
there is no reason to think that the best agreement with the observed profile has to be found 
in the absence of a magnetic field. It is likely enough that the best agreement could be 
obtained in the presence of a deterministic or microturbulent magnetic field, using different 
values of the parameters. For example, in panels c and d of Figure [6] we show the best 
theoretical profiles that we have obtained in the presence of a microturbulent magnetic field 
of 5 G and a vertical magnetic field of 40 G, choosing different values of the free parameters 
(see caption to Fig. [6]). 
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Fig. 6. — Panel a: observed Q/I profile of the Ba II D 2 line ( Stenflo fc Keller 1997h . Panel b: theoretical 
Q/I profile obtained in the absence of a magnetic field, choosing a Doppler width (AAd) of 30 mA, an 
anisotropy factor (w) of 0.037, and a continuum ef/(£j) m ax of 9 x 10~ 5 . Panel c: theoretical Q/I profile 
obtained in the presence of a microturbulent magnetic field of 5 G, choosing the following values of the free 
parameters: AAr>=30 mA, w=0.052, ef/(e l j) max —l.77 x 10~ 4 . Panel d: theoretical Q/I profile obtained 
in the presence of a vertical magnetic field of 40 G, choosing the following values of the free parameters: 
AA D =30 mA, w=0.037, £ / c /(4)ma X = 7.5 x 10~ 5 . 



4.2. D2 Line - The Influence of a Magnetic Field on the Emergent Polarization 

Depending on its strength, and on its direction with respect to the local vertical and 
to the direction of the scattered radiation, a magnetic field will differently modify the linear 
and circular polarization of the line through the Zeeman and the Hanle effects^. It is well 
known that the Zeeman effect produces in general elliptical polarization, which degenerates 
into linear polarization if the magnetic field lies on the plane perpendicular to the line-of- 
sight (LOS), and into circular polarization if the magnetic field lies along the LOS. The 
Zeeman effect dominates the polarization of the scattered radiation if the splitting among 



8 Hereafter, regardless of the particular regime (Zeeman effect regime, incomplete or complete Paschen- 
Back effect regime) , any polarization signal that originates from the splitting among the magnetic sublevels 
will be referred to as Zeeman effect. 
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the magnetic sublevels is of the same order of magnitude or larger than the Doppler width 
of the line. This criterion gives a critical value of the magnetic field strength for the Ba II 
D 2 line of about 3000 G. However, if the magnetic field is not too weak, and if there are 
no other mechanisms that dominate the polarization, it is possible to identify Zeeman effect 
signatures on the fractional polarization profiles even for intensities much smaller than the 
critical value. As we will see below in Figures [7J and [BJ for the line under investigation 
magnetic fields of about 50 G are enough in order the transverse Zeeman effect to produce 
appreciable modifications of the linear polarization signal. 

On the other hand, as described in § [3j a magnetic field is able to modify the atomic 
polarization, and therefore the polarization of the scattered radiation. Depending on the 
configuration of the magnetic field, and on the geometry of the scattering event, different 
signatures of the Hanle effect can be produced on the polarization profiles. As the observed 
scattering polarization in the Ball D 2 line is dominated by the atomic polarization of the 
upper level, recalling that for the 2 P%/2 level tuf e ~ l/A u g m 10~ 8 s and = 1.33, applying 
equation (flOl) . we find that the line is expected to be sensitive to the Hanle effect for magnetic 
field strengths ranging approximately between 1 G and 100 G (values which are smaller than 
the ones needed for the transverse Zeeman effect to be appreciable). 



In this section we consider the effect of a vertical magnetic field on the theoretical Q/I 



profilqj of Figure [6] (panel b). The results are shown in Figure [7J The first interesting 
feature is the enhancement of the linear polarization at the wavelength positions of the two 
dips between the line-core peak and the two secondary peaks, for magnetic fields relatively 
weak (less than 50 G), for which the influence of the transverse Zeeman effect is negligible. 
In order to understand which physical mechanism is at the origin of this and other features 
shown by these Q/I profiles, in the various ranges of magnetic field intensity, we try to 
distinguish which polarization properties of the emergent radiation are due to the atomic 
polarization effects and which ones are due to the Zeeman effect. To this aim, we can obtain 
interesting information by plotting the profiles obtained through equation (fT5|) according to 
two different strategies: 

a) taking the nominal values for the spherical statistical tensors, but setting B = when 
calculating the energy eigenvalues and eigenvectors, 



9 Hereafter by 'theoretical Q/I profile' we will always refer to the profile obtained by applying equa- 
tion (TT5|) . 



4.3. D2 Line - Vertical Magnetic Field 
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Fig. 7. — Theoretical Q/I profile of the Ba n D2 line in the presence of a vertical magnetic field {Ob =0°). 
In the left panel the magnetic field varies between and 100 G while in the right panel it varies between 
and 1000 G. Note the difference in the scale of the two figures. 

b) setting equal to zero all the spherical statistical tensors, except p[j, but taking properly 
into account the influence of B on the energy eigenvalues and eigenvectors. 

In the former case (Fig. [SJ middle column) we are taking into account only the effects due 
to atomic polarization, neglecting the Zeeman effect, while in the latter case (Fig. [HJ right 
column) we are taking into account only this second effect, within the framework of the 
Paschen-Back effect theory. 

From Figure [HJ, it is clear that the above-mentioned enhancement of the linear polar- 
ization at the wavelength positions of the dips, in the presence of a weak vertical magnetic 
field, is not due to the Zeeman effect, but to the Hanle effect acting only on the isotopes 
with HFS (see § 14.11) . We note in fact that the central peak, which is due to the isotopes 
without HFS is not sensitive, in this range, to the magnetic field. Actually this particular be- 
haviour can be explained in terms of two different mechanisms: the anti-level-crossing; effect 
(briefly introduced in § I2.ip . and the change of coupling scheme of the atomic system! 10 !. We 
note first that if the magnetic field lies along the symmetry axis of the radiation field, only 
the statistical tensors p(j and p§ are different from zero (see also § 12. 3p . As the incomplete 
Paschen-Back effect regime is reached, the HFS magnetic sublevels with the same / quantum 
number separate from each other (see § 12.11) . as a consequence the terms Pq(F,F') (which 
quantifies the corresponding quantum interferences) decrease, and this causes an increase of 
the polarization of the scattered radiation (anti- level-crossing-effect, see LL04 for details). 



10 According to our definition, these effects should be better considered as particular cases of the Hanle 
effect. 
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Fig. 8. — Left column: theoretical Q/I profiles of the Ba II D2 line obtained in the presence of a vertical 
magnetic field. Middle column: theoretical Q/I profiles obtained neglecting the Zeeman effect (only atomic 
polarization effects). Right column: theoretical Q/I profiles obtained neglecting atomic polarization effects 
(only Zeeman effect). The zero of the wavelength scales is taken at 4554 A. 



On the other hand, as already stated in § 12.11 going from the Zeeman effect regime to the 
complete Paschen-Back effect regime, the magnetic field produces an energy eigenvectors 
basis transformation. This transformation implies a changing of the coupling scheme of 
the atomic system, which affects the polarization state of the atomic syste n0. Both these 
mechanisms begin to play an appreciable role as the upper level of the D2 enters the incom- 
plete Paschen-Back effect regime. The order of magnitude of the magnetic field strength 
needed to reach this regime can be estimated from the relation 0.1 < AA^/AAhfs ~ 1, where 
A\ B is the splitting induced by the magnetic field, and AAhf s is the wavelength separation 



11 The complex mechanism of inhibition of atomic polarization transfer discussed in 
a particular consequence of this coupling scheme transformation. 
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( 20021) is 



-23- 




Fig. 9. — Theoretical Q/I profiles of the Ba II D2 line in the presence of a horizontal magnetic field, 
perpendicular to the line of sight (#3 =90°, ^ B = ±90°) 



between the HFS F-levels. Applying this relation to the Ba II D2 line, we find that this 
effect is expected to take place for magnetic fields larger than 10 G, as it is observed in 
Figure [HJ The possibility of an enhancement of the scattering polarization in the presence 
of a vertical magnet i c field , through this kind of mechanisms, was already pointed out by 



Trujillo Bueno et al.l (120021 ) for the case of the Nal D2 line. 



Increasing the magnetic field strength, the transverse Zeeman effect eventually becomes 
appreciable and, besides the previous effect, we see an increase of the polarization at the 
wavelength position of the two peaks on the wings of the profile (see Fig. [7J and Fig. [8]). 
Going to magnetic fields of about 200 G or stronger we enter the transverse Zeeman effect 
regime, and the linear polarization profile takes the typical symmetrical shapd^l (see Fig. [7J 
and Fig. ED. 



4.4. D 2 Line - Horizontal Magnetic Field, Perpendicular to the Line of Sight 

In the presence of a weak horizontal magnetic field, perpendicular to the line of sight, as 
the field increases we observe a decrease of the linear polarization at the wavelength position 
of the central peak, due to the Hanle effect. However, consistently with the fact that the 
Hanle effect has to vanish in the far wings of the line, the two peaks on the wings remain 
almost unaffected, as can be seen in the left panel of Figure [91 Similarly to what happens 
in the presence of a vertical magnetic field, going to intensities of about 50 G or stronger we 
enter the transverse Zeeman effect regime, and the Q/I profile takes the well known shape 



12 Note that here we are plotting the ratio eq{X)/ei(X) and not just sq(X). 
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Fig. 10. — Theoretical Q/I (left panel) and U/I (right panel) profiles of the Ba II D2 line in the presence 
of a longitudinal magnetic field (9b =90°, xb = 0°). For xb — 180° the U/I profile would be the same except 
for an overall sign switch. 

shown in the right panel of Figure M 

4.5. D2 Line - Horizontal Magnetic Field, Directed Along the Line of Sight 

In the presence of a longitudinal magnetic field of increasing strength, there is again a 
decrease of the linear polarization at the wavelength position of the central peak, due to the 
Hanle effect, while the two peaks on the wings are not affected. In this geometry the Zeeman 
effect does not modify the linear polarization and, going to stronger magnetic fields, we enter 
a regime of saturation, as shown in the left panel of Figure [TUl Because of the Hanle effect, 
we have in this case a rotation of the plane of linear polarization. This implies the presence 
of the non zero U/ 1 signal shown in the right panel of Figure [10J Finally we have a typical 
antisymmetric V/I signal due to the longitudinal Zeeman effect, as shown in Figure [TTJ For 
weak magnetic fields (of the order of about 100 G or weaker) the signal increases linearly 
with the magnetic field strength. Going to stronger fields the linearity is slowly lost, and the 
profile starts to saturate. 

4.6. D 2 Line - Random- Azimuth Magnetic Field 

In this section we present the results obtained for the fractional polarization in the 
presence of magnetic fields with a given inclination and a random azimuth (i.e. the results 
obtained in the presence of a magnetic field of given strength and inclination, averaged over 
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Fig. 11. — Theoretical V/I profiles of the Ba II D2 line in the presence of a longitudinal magnetic field. 



the azimuth). Figure [T21 shows the theoretical profiles obtained in the presence of random- 
azimuth magnetic fields with inclinations of 30°, 60° and 90°. In the 30° case we observe 
that increasing the magnetic field strength the linear polarization decreases at the wavelength 
position of the central peak and of the two dips close to it, because of the Hanle effect. For 
magnetic fields of about 100 G it is possible to observe the first signatures of the Zeeman 
effect, which dominates the linear polarization as we further increase its intensity. Note 
that as the main component of the magnetic field is vertical, the Zeeman effect produces 
the typical three lobes profiles with the same signs as in the case of a deterministic vertical 
magnetic field. Similar considerations hold for a random-azimuth magnetic field with an 
inclination of 60° and 90°. As the main component of the magnetic field is now horizontal, 
as far as the Zeeman effect starts to dominate the polarization (which happens for fields of 
about 200 G or stronger in the 60° case, and for fields of about 50 G or stronger in the 90° 
case), we obtain the well known Zeeman effect profiles with the same signs as in the case of a 
deterministic horizontal magnetic field, perpendicular to the LOS. Stronger magnetic fields 
are needed for the Zeeman effect to be appreciable in the presence of a random-azimuth 
magnetic field with an inclination of 60° because in this case the vertical and the horizontal 
components of the magnetic field (which produce Zeeman effect profiles with opposite signs) 
are comparable. 
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Fig. 12. — Theoretical Q/I profiles of the Ba n D2 line in the presence of a random-azimuth magnetic field 
for three different inclinations: 30° (top), 60° (middle), and 90° (bottom). 



4.7. D 2 Line - Microturbulent Magnetic Field 

Averaging the emission coefficient over all the possible orientations of the magnetic 
fields, we can investigate the polarization properties of the line in the presence of a unimodal 
microturbulent magnetic field. As shown in Figure [T5| we observe that the linear polarization 



See Appendix [C] for the details implied in performing this average. 
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Fig. 13. — Theoretical Q/I profiles of the Ba II D2 line in the presence of a microturbulent magnetic field. 

at the wavelength position of the central peak decreases, while the two peaks on the wings 
remain constant as the magnetic field strength is increased. This behaviour can be easily 
understood: in the presence of a microturbulent magnetic field there is no observational 
polarization signal due to the Zeeman effect, while the Hanle effect produces a decrease of 
the linear polarization only in the line-core. For B > 100 G, we enter a saturation regime and, 
contrary to the case of a longitudinal field, we still have a non-zero signal in the line-core. 

4.8. The Polarization of the D\ Line 

To understand the physical origin of the polarization signal observed in the Di line and, 
in particular, the role of HFS, we start our investigation considering separately the isotopes 
138 (without HFS) and 137 (with HFS), before taking into account the contribution coming 
from all the seven isotopes together. We use the values of the average number of photons 
and of the anisotropy factor given in § 12.3} and a Doppler width of 30 mA. In Figure Q3] we 
can see that the emission coefficient Eq of the isotopes without HFS is constant in wave- 
length and equal to zero, while it shows an anti-symmetrical profile, that goes rapidly to zero 
moving away from the line-core, in the isotopes with HFS. For these isotopes the frequency 
integrated Stokes Q emission coefficient is equal to Z e r d3 The profile of tie ratio e^/e', of 
the isotopes with HFS, contrary to the case of the D2 line, goes slowly to zero moving away 
from the line-core. For this reason the effect of the continuum is less important in this line, as 
it just 'pushes' more rapidly to zero the wings of this profile without modifying significantly 



All these properties can be derived analytically, and are briefly discussed in Appendix [Bl 
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138 137 All isotopes 




Fig. 14. — Theoretical profiles of the emission coefficients ej(A) and £q(X), and of the ratio £q(X) / (ej(A) + 
ef), of the Ba II Di line, for the isotope 138, for the isotope 137 and for all the seven isotopes together, in 
the absence of a magnetic field. The emission profiles are normalized to the maximum value of the intensity 
emission profile, calculated taking into account the contribution of all the seven isotopes, (£j) max - The 
vertical dashed lines show the positions of the three groups of transitions (see § I4.ip . The last row shows the 
£q(A)/(e|(A) +ef) profiles without continuum (solid), and with a continuum ef/(e l I ) max of 9 x 10^ 5 (dot). 



its shape. In the following description of our investigation on the magnetic sensitivity of this 
line we have neglected the contribution to the intensity coming from the continuum. 

As seen in Figure [IH within the framework of our modeling approach for weak magnetic 
fields (B < 5 G) it is not possible to obtain the symmetric Q(X)/I(X) profile observed by 



Stenflo et al.l (120001 ) . The theoretical profile that we have obtained has no evident symme- 



tries, and its main peak does not coincide in wavelength with the central peak of the observed 
profile. Changing the Doppler width and the anisotropy factor we can modify the width and 
the amplitude of the peaks, adding a continuum contribution to the Q Stokes parameter we 
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Fig. 15. — Theoretical profiles for the ratio £q/£j oftheBallDi line, in the presence of a vertical, horizontal 
perpendicular to the line-of-sight, longitudinal and microturbulent magnetic field of various intensities. 



can shift the profile along the polarization scale in order to get values of Q/I in the wings 
close to the observed ones, but we cannot really modify the shape of the profile. In Figure [TBI 
we show the theoretical s 1 q/s\ profiles in the presence of a vertical, horizontal perpendicular 
to the line-of-sight, longitudinal and microturbulent magnetic field. In complete analogy 
with the D2 line, the profiles are m odified by the combin e d acti on of the Hanle and Zeeman 
effects. As already pointed out by iTrujillo Bueno et al.l (120021 ) and by ICasini et al.l (120021 ) 
for the Nal Dx line case, it appears that, at this level of approximation, the only way to get 
a symmetri c profile, centered a t the wavelength position of the central peak of the observed 
profile (see IStenflo et al.ll2000l ). is to be in the presence of a magnetic field strong enough 
to enter the transverse Zeeman effect. However, as expected, the transverse Zeeman effect 
produces in the £ 1 q/s\ profile two wing lobes that are more significant than the central one, 
which is not the case of the observed profile. 



-30- 



5. CONCLUSIONS 

The most interesting general conclusion of our theoretical investigation on the magnetic 
sensitivity of the D-lines of Ba II is that the observation and modeling of the Hanle and 
Zeeman effects in these resonance lines provide a novel diagnostic tool for mapping the 
magnetic fields of the upper photosphere and lower chromosphere. 

In particular, the Ball D2 line at 4554 A is particularly interesting because the emergent 
linear polarization has contributions from different isotopes, contributions that are easily 
resolved and have a different behavior in the presence of a magnetic field. As a result, there 
is a differential magnetic sensitivity of the emergent linear polarization at line center (where 
the signal is produced by the even isotopes without HFS) with respect to the line wings 
(where the signals are produced by the odd isotopes with HFS). For instance, for the case 
of a vertical magnetic field with a strength between 10 G and 100 G, approximately, only 
the isotopes with HFS are sensitive to the Hanle effect, which produce an enhancement of 
the scattering polarization at the two Q/I wing wavelengths. For the case of a horizontal 
field between about 1 and 100 G the most conspicuous observable effect is the line core 
depolarization produced by the Hanle effect of the barium isotopes devoid of HFS. In both 
cases, the transverse Zeeman effect begins to play an increasingly dominating role for field 
intensities larger than 100 G, approximately. Useful information on the magnetic sensitivity 
of the Q/I profile of the calculated emergent radiation in the D 2 line can be seen in Figure WI\ 
which corresponds to the case of a random azimuth magnetic field with a fixed inclination. 
Of particular interest is the case of an unimodal microturbulent and isotropic magnetic field 
(see Figure [T3|) . for which there is no contribution from the Zeeman effect and Stokes U and 
V are zero. 

Concerning the enigmatic Ba II Dx line it is important to note that in the absence of 
magnetic fields only the 18% isotopes with HFS are capable of producing linear polarization 
through bound-bound transitions. As with the sodium Di line, this is possible thanks to 
the fact that in the absence of depolarizing mechanisms only the upper and lower levels 
of the Dx line transition in the odd isotopes are significantly polarized. Interestingly, a 
Q/I profile with a conspicuous blue-shifted peak is obtained if only the selective emission 
of polarization components that results from the upper-level polarization are taken into 
account (see the last panel of the third row of Fig. [T4"li^l . Under such circumstances one 
could argue that a detailed radiative transfer solution for the Ba II Di line including the 
Doppler shifts caused by the convective motions and waves that are present in the solar 



15 Note, however, that in the absence of a magnetic held the theoretical Q/I max profile of the Ban Dx line 
(see the last panel in the second row of Fig. [TH) has an antisymmetrical shape, as already found for the Nal 
Di line (see Fig. 2 in lTruiillo Bueno et al.ll2002l ), when only selective emission processes are considered. 



31 



atmospheric plasma could perhaps produce a symmetric Q/I profile for the Ball Di line, as 
observed by Stenflo et al. (2000). However, as shown in this paper, for the Ball D\ line we 
should expect also a significant contribution from "zero-field" dichroism -that is, from the 
selective absorption of polarization components that results from the lower-level polarization. 
In fact, when both selective emission and absorption processes are taken into account through 
the approximation of equation (fl3|) . we then obtain the nearly antisymmetric Q/I profile 
of Figure [THJ Note that there is no possibility of destroying the lower level polarization 
without sim ultaneously destroying the atomic polarization of the upper level of the D± line 
(see Fig. 1 of lTrujillo Bueno et al.ll2002l . and our Fig. [3]for barium levels). As far as dichroism 
is neglected, within the framework of our present modeling assumptions one might then be 
tempted to conclude that the only possibility of obtaining a symmetric Q/I peak for the Ba II 
~Di line is via the transverse Zeeman effect, even though, through this kind of mechanism, a 
profile with wing lobes more significant then the central one is obtained, and magnetic fields 
quite intense are needed. Our approach neglects, however, the radiative transfer effects that 
we certainly have in the real solar atmosphere, with its vertical stratification, horizontal 
inhomogeneities and the Doppler shifts caused by the above-mentioned upflows, downflows 
and waves. Therefore, detailed radiative transfer simulations using realistic solar atmospheric 
models are urgently needed in order to be able to conclude whether a symmetric Q/I profile 
for the Ball Di line with a significant line-center peak may be obtained within the framework 
of the density matrix theory we have applied in this paper, either because of the influence 
of the atomic level polarization of the 18% of the barium isotopes endowed of HFS (which 
would require the presence of atmospheric regions with very weak fields), or due to the 
transverse Zeeman effect of all the barium isotopes (which would require the presence of 
a sufficiently strong magnetic field, and the absence of significant saturation effects at the 
line-center wavelength) . 

Finally, we would like to finish this paper by emphasizing the importance of pursuing 
high-spatial resolution polarimetric observations of the Ball D 2 line (e.g., via Fabry- Perot 
polarimetry) in order to help decipher the spatial and temporal fluctuations of the magnetic 
field vector, in both active and quiet regions of the solar atmosphere. 
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Fig. 16. — Le/t column: theoretical profiles for the ratios (from the top) £q/£j, T]q/t]j and for their 
difference, in the presence of a vertical magnetic field of various intensities. Right column: line profiles taking 
into account the continuum contribution to the intensity emission and absorption (£//(4)max = 9 x 10 -5 , 
■qf calculated through eq. |A2| ). Calculations refer to the Ba II D2 line. 



A. RESULTS OBTAINED TAKING INTO ACCOUNT DICHROISM 

EFFECTS 

We show here the results obtained taking into account the absorption effects in the D 2 
line. Including the contribution of the continuum, equation (|13p takes the form 

„ eirfo") _ qjtfo") (A1) 
~ ej(i/,fi) + ef ri^M+v! ' 

where 77^ is the continuum contribution to the total intensity absorption. This quantity can 
be calculated from the continuum intensity emission coefficient through the relation 

Vt = ; — r , (A2) 

11 Bp(v ,T = 5800K) ' v ; 
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Fig. 17. — Same as Figure [TH] for larger values of the magnetic field. 



where Bp(vq,T = 5800K) is the Planck function calculated at the central wavelength (i/q) of 
the line we are investigating and at the effective solar temperature. Figures [TBI and fTTl show 
the results for the D 2 line. We consider first the results obtained without continuum. In the 
absence of magnetic fields, since the atomic polarization in the lower level is much smaller 
than in the upper level of D 2 (see Fig. EJ), the absorption effects are completely negligible, as 
already observed in §|U Introducing a vertical magnetic field the profile of the ratio i]qIt]\ is 
modified by the transverse Zeeman effect and, as the magnetic field is increased, it assumes 
a shape which is very similar to the one observed for Sg/ej (see Fig. [IT]) . For this reason, 
the profile obtained through equation ( {TBI does not show any detail due to the transverse 
Zeeman effect since the contribution coming from emission and absorption cancel out. The 
situation is somewhat different in the presence of the continuum. In this shown in 

Figure [TTJ for magnetic fields stronger than about 100 G the contribution coming from the 
absorption term becomes more important. At this point it is important to remember that 
the continuum has been considered as a parameter in our investigation and, as stressed in 
§ 14.11 we have to be careful when dealing with spectral details that find their origin in this 
physical aspect of the problem. Similar considerations can be done for the other magnetic 



4934 4934.2 

A (A) 

Fig. 18. — Theoretical Q/I profile of the Ba II Di line calculated according to equation [T3l in the absence 
of magnetic fields and without any contribution of the continuum. For the Doppler width, the anisotropy 
factor, and the average number of photons have been used the same values of § 14.81 

field geometries considered in this work. 

Figure [TBI shows the ratio (sq/ej) — (t]q/t)][) for the Di line in the absence of magnetic 
fields. As expected, since the atomic polarization is quite similar in the upper and lower 
levels of the Di, dichroism is much more important in this line, and the profiles obtained 
through equations (fT4"|) and ( fi~3l) are quite different from each other. However, even taking 
into account dichroism, we are not able to reproduce the observed profile. 

B. ANALYTICAL RESULTS FOR A TWO-LEVEL ATOM APPLIED TO 

THE Bail D LINES 

In this work we have described the Ba II ion through a three level model atom that 
allowed us to study both the Dx and D 2 lines; we have taken into account the ground level 
polarization, the stimulated emission effects and, finally, the effects of a magnetic field. At 
this level a numerical approach of the problem is absolutely necessary. However, as described 
in detail in § 10 of LL04, by introducing some simplifying approximations it is possible to 
obtain analytical expressions for the atomic density-matrix elements and, through these, for 
the radiation transfer coefficients of the atomic system. These analytical expressions are very 
useful in order to understand the physics of the phenomenon under investigation, which could 
remain quite hidden within a numerical approach. The basic approximation is to consider 
a two-level atom. The only difference between our three-level model atom and a two-level 
atom lies in the fact that in the SEEs of our model the ground level feels the effect of both 
transitions towards the two upper levels considered. However, the numerical solution of the 
SEEs showed that the upper level of the D 2 line is much more polarized than the ground 



-35 - 



level, so that, as far as the D2 line is considered, the polarization of the ground level, as a 
first approximation, can be neglected and our atomic system can be treated as a two-level 
atom. We can therefore apply the equations of § 10 of LL04 to calculate the upper level 
density-matrix elements or the emission coefficients of the various Ball isotopes, described 
as two-level atoms with HFS, and two-level atoms without HFS, with unpolarized lower 
level. Neglecting the stimulated emission effects (which is a good approximation whenever 
the incident radiation field, as in our case, is weak), and in the absence of the magnetic field, 
the emission coefficient of a two-level atom without HFS is (see eq.[10.16] of LL04) 

ei(u, Q) = ^-AfeB(a e Je -> a u J u )(f)(is - v) 



x 



J2WK(Je,Ju)(-l) Q T«(i,n)J« Q (v ) , (Bl) 



KQ 



where Aft is the number density of atoms in the ground leve0, and where 



1 IK" 2 



W K (J e ,J u ) = 3(2 J u + l)\ T T T \ ■ (B2) 

For 90° scattering of a radiation field with cylindrical symmetry around the direction of 
propagation (as in our case), from equation fIBip we obtain 

= e Q 3^(1/2,3/2) 
PQ ~ sj -4/ W -W 2 (l/2,3/2) ' {B6) 

This expression shows that for the isotopes without HFS the fractional polarization pq does 
not depend on frequency, as found in § 14.11 (Fig. EJ). Substituting the numerical values of the 
various quantities (w = 0.037 and ^(1/2, 3/2) = 0.5) we obtain the value of 1.4%, as found 
with numerical calculations. The expression of the emission coefficient of a two-level atom 
with HFS (eq. [10.166] of LL04) is much more complicated and will not be written here. 
Anyway, as shown in § 10.22 of LL04, at frequencies very distant from the 'center of grav- 
ity', the multiplet behaves in resonance scattering as a simple transition between two levels 
without HFS. The fractional polarization, therefore, at these frequencies is still described by 
equation (1B3I) . This result justifies the fact that at large distances from the line center the 
fractional polarization of the isotopes with HFS reaches the same value as the isotopes with- 
out HFS. It is important to stress that this asymptotic behaviour of the isotopes with HFS is 
strongly dependent on the interferences between the various HFS magnetic sublevels of the 



16 Note that for a two level atom M = M e +M u while, for our model atom, TV = Me + M® 2 +M^ >1 . As far 
as lower level polarization is neglected, this is the only difference between a two-level atom and our model 
atom. 
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D2 upper level. Because of the small frequency distance between the various components of 
the HFS multiplet with respect to their natural width, the analytical expression of the emis- 
sion coefficients of the isotopes with HFS cannot be simplified in the neighbourhood of the 
various transitions. Nevertheless, we can qualitatively justify the decrease of the fractional 
polarization at the wavelength positions of the various components of the HFS multiplet by 
considering the frequency integrated emission coefficients of the isotopes with HFS 



JAu 



71 KQ 

where the interval Av is sufficiently broad to fully cover all the Zeeman components of the 
line, and where 

[W K ( y aiJiI,a u J u )}Ms = W K (Je,Ju)[DK((y u J u I)]Ms ■ (B5) 
The quantity [D ic(oi u Ju^)]ih is t ne depolarizing factor due to HFS and it is given by 

[DK(aJI)U = (^£(2F+1)(2F' + 1){^ J p , f } 2 

x l —— . (B6) 

1 + 2mis a j IF/ia j IF /A(aJ -> a t J t ) 

For the isotope 137 the depolarizing factor [D 2 (J=3/2, J=3/2)] hfs is equal to 0.27, and the 
frequency integrated fractional polarization, that in complete analogy with equation ( 1B3I ) 
has the form 

~ _ £q _ 3W 2 (Je, Ju)[DK{u u JuI)\hfs 

PQ ~ ej ~ A/w - W 2 (J e , J u )[D K (a u J u I)} hfs ' 1 ' 

is equal to 0.0038. Comparing this value with 0.014, the value of pq previously found for 
the isotopes without HFS (note that for these isotopes pq is equal to Pq), we can clearly see 
the depolarizing effect of the HFS. 

For the Di line the quantity W^(l/2, 1/2) is zero and, taking into account equations (IBID 
and (]B7p . it is easy to see that the Stokes Q emission coefficient of the isotopes without HFS 
has to be constant and equal to zero, while the frequency integrated Stokes Q emission 
coefficient of the isotopes with HFS has to be zero, as found with numerical calculations 
in § 14.81 However, as previously said, for the case of the Di line the approximation of 
unpolarized lower level is not good anymore. Some analytical results that it is possible to 
obtain taking into account the polarization of the lower level are derived in LL04. 

In the presence of a vertical magnetic field equation (1B7I) generalizes into 

PQ ~ ij ~ A/w - [W 220 (a e J £ Ia u J u ; B)] Ms ' 1 j 
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B(C) 

Fig. 19. — Frequency integrated fractional polarization for the isotope 137 in a 90° scattering event in the 
presence of a vertical magnetic field, as function of the magnetic field strength. Calculations refer to the 
Ba II D 2 line. 



where the general expression of the quantity \WKK'Q{citJfIa u J u ] B)}^ is given in LL04 
(eq. [10.167]). Figure [191 shows that there is an increase of Pq for vertical magnetic fields 
ranging between and 500 G (see Fig. [7]). 



C. AVERAGE OF THE EMISSION COEFFICIENTS OVER THE 

MAGNETIC FIELD DIRECTIONS 



The only quantity that depends on the magnetic field orientation in the general expres- 
sion of the emission coefficients (§ I2.1[ eq. [6]) is the product 



B 



auJuI Po:(K,K) 



B 



(CI) 



where the label B means that the quantity is calculated in the magnetic field reference system. 
We already observed (§ [3]) that the spherical statistical tensors calculated in the magnetic 
reference system, because of the symmetry of the problem, do not depend on the azimuth of 
the magnetic field, xb, but only on its inclination with respect to the local vertical, 9b- For 
this reason all the dependence on xb is included into the geometrical tensor [TJf(j, 
The expression of the tensor Tq in terms of rotation matrices is given by (eq. [5.159] of 
LL04) 

T Q^ 0)1 R = J2 tK P^ V URo) , (C2) 

J B 

p 

where T> is the rotation matrix, Rq is the rotation bringing the reference system (e (f2), e&(Q), fi) 
with e a (fi) the reference direction, into the reference system with the z axis directed along 
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the magnetic field, and where tp(j) is a scalar quantity that does not depend on the partic- 
ular geometry of the problem (cf. eq. [5.160] of LL04). The relation between the geometrical 
tensor Tq calculated in the magnetic field reference system, and the same quantity calcu- 
lated in the local vertical reference system (with the z axis directed along the local vertical) 



is 



E 



rum 



V 



AQ 



(Rl 



(C3) 



where the label V means that the corresponding quantity is calculated into the local vertical 
reference system, and where Ri is the rotation bringing the local vertical reference system 
into the magnetic field reference system. In the geometry of our problem the rotation Ri 
is defined by the Euler angles (xb,#b,0). Obviously, because of its definition, [Tj^(j, ft)] v 
does not depend on the magnetic field orientation: in analogy with equation flU2l) we can 
write 

T A K (m} v = j2tm^p A (R2) , (C4) 



where R<i is the rotation that brings the reference system (e a (fi), e^fi), f2) into the local 
vertical reference system. Substituting equation (1C3I) into the expression (IClj) we obtain 



1 Q 



ti, n) 



J B L 



I Jul sJ^'< 



B 



E 

A 



B 



(C5) 



where we see that all the dependence on xb is included in the rotation matrix V^q(Ri). 
Recalling the expression of the rotation matrices in terms of the reduced rotation matrices 
(cf. eq. [2.68] of LL04), and recalling the Euler angles corresponding to the rotation Ri, we 
have 

V« Q {R 1 )=e-^d K AQ {6 B ) , (C6) 

where d A g is the reduced rotation matrix, which depends only on the second Euler angle 
of the rotation, in our case the inclination of the magnetic field. Averaging on xb reduces 
therefore to calculate the integral 

"2- 



1 

2^ 



-iAxi 



dXB 



(C7) 



It is easy to see that the integral is different from zero only if A — 0, and that, in this case, 
it is equal to 1. Averaging the emission coefficients over xb, the magnetic field azimuth, is 



therefore equivalent to substitute expression (IClj) with 



van) 



J B 



Now all the quantities depend only on the magnetic field inclination, and we can complete 
numerically (for example by a Gaussian quadrature) the average over the magnetic field 
orientation. 
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